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Abstract 

We present a systematic semiclassical procedure to compute the par- 
tition function for scalar field theories at finite temperature. The central 
objects in our scheme are the solutions of the classical equations of mo- 
tion in imaginary time, with spatially independent boundary conditions. 
Field fluctuations - both field deviations around these classical solutions, 
and fluctuations of the boundary value of the fields - are resummed in a 
Gaussian approximation. In our final expression for the partition func- 
tion, this resummation is reduced to solving certain ordinary differential 
equations. Moreover, we show that it is renormalizable with the usual 
1-loop counterterms. 
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1 Introduction 

Finite-temperature field theory [1] is the natural framework for the study of 
phase transitions, and of the thermodynamic properties of equilibrium states. 
Applications range from the investigation of the phase structure of the strong 
and electroweak interactions, and the related applications to the early universe, 
to the low-energy effective field theories in particle physics and condensed matter 
systems. 

However, finite-temperature field theories often face a major difficulty: the 
plain perturbation expansion [5-8] is ill-defined due to the presence of infrared 
divergences in the bosonic sector, and often gives meaningless results. In the 



case of hot QCD, for instance, one can say that the domain of validity of the 
naive weak-coupling expansion is the empty set [2]. This challenge stimulated 
the development of resummation techniques that reorganize the perturbative se- 
ries, and resum certain classes of diagrams, thereby improving the perturbative 
expansion (see [3,4] for recent reviews). Some of these techniques amount to 
using an effective theory in order to separate the scales T, gT, and g 2 T [9-14]. 
Others use modified quasi-particles as the starting point of the perturbative 
expansion [15], leading to a significant improvement of the convergence of the 
expansion when the mass of these quasi-particles is properly chosen (one can 
also mention Rcfs. [27,28], where a simple phcnomcnological model of massive 
quasi-particles was successfully used in order to reproduce the pressure of the 
quark-gluon plasma obtained in lattice simulations). In other approaches that 
aim at maintaining thcrmodynamical consistency, one reorganizes the perturba- 
tive expansion of the thermodynamical potential around two-particle irreducible 
skeleton diagrams [16-21]. Finally, some of these techniques are based on a sys- 
tematic use of the Hard Thermal Loop effective action [22-26], i.e., on the 
assumption that the HTLs provide a good description of the quasi-particles in 
the plasma, and of their interactions. 

A somewhat different approach, which can also be interpreted as a resumma- 
tion of an infinite set of perturbative diagrams, is provided by the semiclassical 
approximation [29,30]. Since the partition function of a given system can be 
cast in the form of a path integral whose weight is the exponential of minus the 
Euclidean action, an expansion around Euclidean classical solutions is quite nat- 
ural. This program has been carried out in the case of one-dimensional quantum 
statistical mechanics for particles in a single-well potential in [31], and also for 
double wells [32] . From the mere knowledge of the classical Euclidean solutions 
of the equation of motion, the full semiclassical series for the partition function 
was constructed 1 . Later, these results were generalized to the case of a particle 
in a central potential in an arbitrary number of dimensions [35]. In both cases, 
excellent results were obtained, for instance, for the ground state energy and 
the specific heat, in the case of the quartic potential. 

In this paper we develop a similar semiclassical procedure to compute the 
partition function for thermal scalar field theories with a single-well potential. 
We expand around Euclidean classical fields, whose value on the boundaries of 
the time interval are taken to be independent of space. These solutions are as- 
sumed to be known to all orders in the interaction potential (either analytically 
or numerically). Then, we incorporate fluctuations around these classical tra- 
jectories, as well as space-dependent fluctuations of the boundary value of the 
classical fields. All these fluctuations are kept only in a Gaussian approximation, 
although it is in principle possible to go systematically beyond this approxima- 
tion. We also provide a diagrammatic interpretation of our results, connecting 
our formalism to the ordinary perturbative expansion, and identifying the classes 

1 The equivalent problem in quantum mechanics at zero temperature was previously studied 
by DcWitt-Morettc [33] for arbitrary potentials, and by Mizrahi [34] for the single-well quartic 
anharmonic oscillator, using similar techniques. For a more complete list of references on the 
semiclassical series in quantum mechanics, see [31] 
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of diagrams that are resummed in our approach. The implementation of the 
renormalization procedure in this semiclassical treatment is discussed in detail 
at the end. 

Since the procedure we propose is infrared finite, we believe it represents an 
interesting alternative to other rearrangements of perturbation theory at finite 
temperature. In this paper, we present the general semiclassical framework 
for an arbitrary potential. We leave the application to the case of a scalar 
theory with quartic self-interactions, and comparisons with other methods to a 
following publication [38]. 

The paper is organized as follows: 

In section 2, we discuss the computation of the partition function, starting 
from its expression in terms of a functional integral, within the semiclassical 
approximation. Although the main result is, of course, well known, we focus 
our discussion on the role played by the boundary conditions in Euclidean time 
at finite temperature. We also recall the diagrammatic interpretation of the 
integration over quadratic fluctuations around the classical field. 

In section 3, we present a systematic procedure to incorporate effects from 
fluctuations of the boundary value of the field in the computation of the partition 
function. We explain how one can perform an expansion in those fluctuations in 
a consistent way, provided that one knows the classical solutions for the problem 
with constant boundary conditions. We derive formulas that incorporate the 
effects of these fluctuations up to quadratic order. These formulas depend only 
on the classical field itself, and on a basis of solutions for the equation of motion 
for small fluctuations around the classical field. 

In section 4, we expand the classical action to second order in the boundary 
fluctuations, and discuss diagrammatically the meaning of this expansion in 
terms of the boundary value of the field. This leads to our final expression 
for the partition function in terms of quantities that can be straightforwardly 
obtained in explicit form for a given potential once one knows the classical 
solution mentioned previously (at least numerically). This expression, however, 
still needs to be renormalized. 

The renormalization procedure, which resembles the usual perturbative pro- 
cedure, is discussed in section 5. There, we show how to obtain a finite expres- 
sion for the partition function through the introduction of only two counterterms 
in the action, plus the subtraction of the zero point energy. 

We present our conclusions and outlook in section 6. Finally, in the ap- 
pendix, we illustrate the procedure in the case of the free theory. As mentioned 
above, the non-trivial example of the quartic potential will be addressed in detail 
in another publication. 

2 Small fluctuations around a classical solution 

We want to calculate the partition function Z = Tre _/3H for a system of inter- 
acting scalar fields, making use of a semiclassical approximation. Our starting 
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point is the expression of Z in terms of path integrals : 



[D<p{x)\ J [D<f>(r,x)} e 

4,(-f3/2, X )=<t,(f3/2, X )=<p{x) 

where S E [(f)] is the Euclidean action of the field: 

+P/2 

S E [<f>]= J drd 3 x 

-f}/2 



(1) 



^d^cp + \m 2 <j> 2 



u(</>) 



(2) 



Assume, for the time being, that we know the solution 4> c (t,x) of the classical 
equation of motion, that takes the value f(x) on the boundaries of the time 
interval : 



(D E + m 2 )0 c (r, x) + U' (Mr, as)) = 



(3) 



where we denote by O e = — (d 2 + V 2 ) the Euclidean D'Alembertian operator. 

A classical solution is a (local) minimum of Se- Next, in the functional 
integration over <((t,x) in eq. (1), one assumes that the integral is dominated 
by field configurations in the vicinity of that classical solution, i.e., by small 
fluctuations around this classical solution. In order to evaluate the integral in 
this approximation, one writes : 

<j>{r,x) = 4> c (r,x) +v(t,x) , (4) 

and expands the Euclidean action to second order in the fluctuation t](t, x) : 

sm = sm 



+\j{d* Xl ) E {t 



x 2 ) 



s 2 s^ 



5<j)(xi)5<j)(x2) 



vfaWxi) + 0(r, 4 ). (5) 



In this equation, we have used the shorthands x = (t,x) and J(d 4 x) E = 

u 

5 2 SM 



$-13/ 2 ^ T i d 3 x. For the sake of brevity, let us also introduce the following no- 
tation : 



r, 1 A^jr, = /(d^ij^a*), 



5<t>(xi)6</)(x2) 



r](xi)r]{x 2 ) , (6) 



4>=4>c 



where A[0 C ] is a symmetric "matrix" that depends on the classical solution </> c 
(with continuous indices spanning [—(3/2,(3/2) x R 3 ). 

The Gaussian functional integration over the fluctuation r\ must be per- 
formed with the constraint that the fluctuation r\(j, x) vanishes at the time 
boundaries, 

Vx, r ] {-P/2,x)=r ] {(3/2,x) = 0, (7) 
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because the classical background field already saturates the boundary condi- 
tions. Let us therefore call A* [cf) c ] the restriction of the operator A[cj) c ] to the 
subspace of fluctuations r\ that obey these boundary conditions. We can write : 



Z « / [D<p{x)\ e 




s ^ [det (A*[<j> c })] 



-1/2 



(8) 



In order to compute the scmiclassical calculation of Z, one must now inte- 
grate over the boundary value of the field, <p{x). However, before we pursue this 
calculation, it is useful to recall the nature of the diagrams that are contained 
in the square root of the functional determinant. It is well known that the 
Gaussian integration over fluctuations above a given background field amounts 
to calculating the one-loop correction to the effective action. However, for this 
correspondence to be valid, one must integrate over all the periodic fields rj(x). 
In our case, the Gaussian integration involves only fields r) that vanish on the 
time boundaries (see eq. (7)), i.e., only a subset of all the periodic fields. There- 

fore, the quantity [det (A*[</> c ])] is a part of the one-loop effective action, 
but does not contain all the terms that would normally enter in the effective 
action at this order 2 . With this caveat in mind, a typical diagram included in 
this quantity is displayed in figure 1, in the case of a field theory with a quar- 
tic coupling. It is important to remember that the propagator represented by 
the dashed line differs from the complete time ordered propagator, because it 
corresponds to a subset of all the periodic modes. 



Figure 1: Typical 1-loop diagram included in the integration over fluctuations 
around the classical solution in the Gaussian approximation. The lines termi- 
nated by a cross denote the classical solution with a fixed boundary condition 
ip(x). The dashed line can be seen as the propagator around the classical field, 
for a fluctuation that vanishes at the time boundaries. 

Both the classical action and the determinant in eq. (8) depend on the field 
<p(x) on the boundary, through the dependence of the classical solution cf> c on 
the boundary conditions in eq. (3). In fact, the classical solution (f> c can be 
represented diagrammatically as the sum of all the tree diagrams terminated 

2 This distinction can also be seen by studying the eigenfunctions of the operator A[</> c ], 
on the space of periodic functions and on the space of functions that vanish at r = ±/3/2 
respectively. 
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by the boundary field <p(x). The easiest way to see this is to write Green's 
formula for the solution of eq. (3). Let us first introduce a Green's function of 
the operator □ + m 2 : 



[8 2 T , + V 2 , - to 2 ] G°(t, x; t', x') = S(t - t')6(x - x') . 



(9) 



This Green's function is not unique, but we can postpone its choice for later. 
Let us multiply this equation by the classical field </> c (r', x'), and integrate over 
t' and x'. This gives : 



/•/J/2 r 

(r, x) = dr' d 3 x' Mr', x 1 ) [0 2 T , + V 2 , - to 2 ] G°(t, x; t', x') . (10) 

J-0/2 J 



Now, multiply the equation of motion for <P c (t',x') by the Green's function 
G°(t,x;t' ,x'), integrate over r', and subtract the resulting equation from the 
previous one. This leads to : 



/pi* r 
dr' / d 3 x' G°{t,x-t',x') U'{4> c {t',x')) 

1-0/2 f _ 

+ / dr' d 3 x / G°(t,x;t / ,x') 

J-0/2 J 



di, + V 2 - di, 



,(r',x'), 
(11) 



where the arrows on the differential operators on the second line indicate on 
which side they act. The second line can be rewritten as a boundary term, by 
noting that : 



A 



B = d^\A 



Of* - df, 



D 



(12) 



In eq. (11), the boundary in the spatial directions does not contribute to the 
classical field at the point x because the free propagator decreases fast enough 
when the spatial separation increases. Thus, we are left with only a contribution 
from the boundaries in time. At this point, since the boundary conditions for 
4> c consist in specifying the value of the field at r' = ±(3/2, while its first time 
derivative is not constrained, it is very natural to choose a Green's function G° 
that obeys the following conditions 3 : 



G°(t, x; -(3/2, x') = G°(t, x; +(3/2, x') = 



(13) 



3 It is in general always possible to impose two conditions on a Green's function of O e +m 2 , 
because the zero modes of this operator form a linear space of dimension 2. The conditions 
of eq. (13) are explicitly realized by : 

G o (Ti x; T , :X , ) = J J0_ e ». ( „-„0 [ e{T _ T>) **Wr - I)) sinh K (.' + I)) 



where we denote cjfc = \/ k 2 + m 2 . 



+ 9(t' - r) 



uj k sinh(aj fe /3) 
ui k sinh(oj fe /3) 
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With this choice of the propagator, we obtain the following formula for <p c (r, x) 



f/3/2 f 

(t,x) = dr' d 3 x / G°(t,x;t',x') U' (4> c {t' ,x')) 

J-0/2 J 



J d 3 x' <p(x') [dr>G a (T,x;T',x r ) 



r'=+/3/2 

T> = -0/2 



(14) 



This formula tells us how the classical solution <fi c depends on the boundary 
value (p(x). If the first term in the right hand side - involving the derivative 
U' of the potential - were not there, then the relationship between 4> c and 
the boundary value ip would be linear. This only happens in a free theory. 
When there are interactions, one can solve eq. (14) iteratively in powers of U'. 
This expansion can be represented diagrammatically by the sum of the tree 
diagrams whose "leaves" are made of the boundary field f(x). An example of 
such a diagram is illustrated in figure 2, in the case of a </> 4 interaction of the 
fields. Notice that, when the boundary field is small 4 , this sum of trees can be 




Figure 2: Diagrammatic expansion of the classical field in terms of the boundary 
value of the field (black dots). 

approximated by the zeroth order in the expansion in powers of U', which is 
independent of the interactions. On the other hand, for large values of tp, it is 
important to keep the full sum of tree diagrams that are summed in <p c , because 
all the terms in the expansion can be equally important. Therefore, we already 
see an important feature of our approximation scheme: although the quantum 
fluctuations are only included at the 1-loop level, it treats the boundary field to 
all orders, allowing a correct treatment even for non-perturbatively large values 
of <p(x). 

By this, we mean that the interaction term is smaller than the kinetic term in the action. 
This condition depends on the particular momentum modes one is interested in. 
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3 Expansion in fluctuations of the boundary 



3.1 Preliminary discussion 

The integration over boundary configurations that remains to be performed in 
eq. (8) makes the semiclassical approximation for Z rather involved; first, we 
must solve the partial differential equation (3) for an arbitrary <p(x), and this 
will not be feasible analytically in general. Even numerically, this is a very 
complicated task. Besides that, the only functional integral over ip that one is 
able to perform analytically is a Gaussian integral. In order to circumvent these 
problems, we are forced to perform some further approximations. 

One can see readily in cqs. (3) that the classical equation of motion reduces 
to an ordinary differential equation in the case where the field ip(x) on the 
boundary is a constant ipo. In this case, the classical solution <f> x (r,x) becomes 
a function 0o (r) of the time only : 

(-d 2 T +m 2 )Mr) + U'(Mr))=0, 

M-0/2) = <M/V2) = • (15) 

Such a simplification of the classical equation of motion would make the problem 
much more tractable by analytical or numerical methods. This remark suggests 
that we decompose the boundary field ip{x) into a constant part ipo, and a 
fluctuation : 

<p(x) = <p + £(x) . (16) 

The solution of the classical equation of motion can therefore be expanded in a 
similar manner : 

<P c (t,x) = 0o (t) + 0i (r, a;) + fo{r,x) H , (17) 

where n is of order n in £ (there are terms of arbitrarily high order in £ if the 
equation of motion is non linear). Having done this, we can rewrite the path 
integral over ip{x) in eq. (8) as follows : 

Za jchpo J [D£{x)\ e- 8 *^ [det (A*[0 C ])]- 1/2 . (18) 
-°° (««))=o 

Notice that the integration over £(x) must be performed with the constraint 
that 

(£(x)) = J d 3 x^x) = 0, (19) 

since the "uniform component" (i.e., the zero mode) of the boundary condition 
is already included in tpo. An unrestricted integration of ^(a;) would therefore 
overcount the contribution of this zero mode. 

In the following, we are going to assume that the first term 0o(t) can be 
determined with an arbitrary accuracy - it can be determined analytically in 
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certain cases, while in general it is obtained by solving numerically an ordinary 
differential equation. Moreover, the dependence on <po will always be treated 
exactly. Only the terms that are of higher order in the fluctuation £ of the 
boundary field will be treated in some approximate way. Doing this allows 
us to preserve the benefits of treating correctly the interaction term when the 
boundary field is large, since only the fluctuations of the boundary field are 
assumed to be perturbative. 

A natural approximation to obtain the dependence on £ is to do a Gaussian 
approximation around £ = 0. As we shall see shortly, in order to find the 
classical action S E [<j> c ] at order two in the fluctuation £ (a;) of the boundary, it 
is enough to obtain the classical solution <f> c at order one in £(x). 

Moreover, to be consistent with the Gaussian approximation for S B [<p c ], we 
only need to evaluate the determinant at lowest order in £(a:), i.e., at order zero. 
Indeed, the Gaussian integration over the fluctuations £ corresponds to a one- 
loop correction in the background (f>o. However, as we have seen in the previous 
section, the functional determinant in eq. (18) is already a one-loop correction. 
Therefore, keeping the £ dependence in this determinant would give higher loop 
corrections when we integrate over £, but only a certain subset of all the 2-loop 
corrections would be included. Doing so is not forbidden by any fundamental 
principle, but it would arguably make the calculation more complicated; and 
moreover this would alter the renormalization of the final result. Indeed, as we 
shall sec later, by not expanding in £ the functional determinant in eq. (18), we 
will eventually obtain an expression whose ultraviolet divergences are precisely 
those of the one-loop effective action. For these reasons, we are going to evaluate 

Zs* J *po [det (A*k* ])r 1/2 / [D&x)] e~ S ^ ] ■ (20) 

(««))=o 

3.2 Correction to <p c due to boundary fluctuations 

The next step is to find the correction <pi(r,x) to the classical solution <fi c . In 
order to find the equation obeyed by <p\, simply replace <j> c by <j>o + 4>i in eq. (3). 
By dropping all the terms that are of order higher than unity in (f>i (since they 
are at least of order two in £), and using the equation obeyed by <j>o, we obtain 
the following (linearized) equation for <f>\ : 

"(□ E +m 2 ) + C7" (Mr)) 

with the boundary condition: 

0! (-0/2, a:) =<M/3/2,x) • (22) 

In the following, we also need the Green's formula for the variation <f>i of the 
classical field. The derivation is very similar to the derivation of eq. (14), and 



= 



(21) 
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we shall not reproduce it here. The main difference compared to eq. (14) is that 
we need a Green's function for the operator (O e + m 2 + U" (((o(r))), 

[d 2 T , + V 2 , - m 2 - U" (Mr'))} G(t, x; t', x') = S(r - r')S(x - x') , (23) 

instead of the free propagator G° that we have introduced earlier. Again, this 
propagator must obey the boundary condition 

G(t, x- -(3/2, x 1 ) = G(t, x- (3/2, x') = 0. (24) 

In terms of the fluctuation £ and of the propagator, the first order correction to 
the classical solution reads: 

Mr,x) = J d 3 x' Z{x')[d T ,G{T,xy,x?)tiZ±%l ■ (25) 

Notice that, since the background field 4>o does not depend on space, the prop- 
agator G depends only on the difference x — x' . Thus, we can get rid of the 
spatial convolution by going to Fourier space: 

Mr, k) - C(k) [d T >G(r, r' , k)\ T ^_% , (26) 
where the propagator in Fourier space is defined by 

[d 2 T , - (k 2 + m 2 ) - U" (M^))] G(t, t', k) = S(r - r') , (27) 

and 

G(T,-(3/2,k)=G(r,P/2,k) = 0. (28) 
3.3 Propagator in the background O 

It is fairly easy to determine the propagator G that obeys eqs. (27) and (28) in 
terms of two linearly independent solutions of the homogeneous linear differen- 
tial equation: 

[8 2 - (m 2 + k 2 ) - U" (Mr))} V(r, k 2 ) = . (29) 

Let r]a(r,k) and rjb(T,k) be two such independent solutions 5 of (29). In order 
to construct from r\ a ^ a solution of eqs. (27) and (28), let us first introduce the 
following object: 

f2(r, r', k 2 ) ee Va (r, k 2 ) Vb (r' , k 2 ) - Vb (r, k 2 )^(r' \ k 2 ) . (31) 
5 When k = 0, it is straightforward to verify that : 



Va(r;0) = 4>o{t) , 




obey (29). However, this construction fails when fc ^ 0. 
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It is trivial to check that 0(r, r', A; 2 ) satisfies eq. (29), both with respect to the 
variable r and to the variable r'. Let us then consider the following quantity: 



H{ , k 2, = n(P/2,r,k 2 )n(T',-p/2,k 2 ) 

fi(/3/2,-/3/2,fc 2 ) 
n(f3/2,T>,k 2 )n(T,-p/2,k 2 ) 
Q((3/2,~p/2,k 2 ) 



H(r,r',k 2 ) 



if t > t' , 



if r < r' 



(32) 



This quantity obeys eq. (27) if r ^ r'. Moreover, although H(t,t' ,k 2 ) is 



continuous at r = r', its first time derivative is not, and one has: 

Jim (0 T ,tf (t, r', fe 2 ) | r , =T+e - a T , H(r, r', fc 2 )| T , =r _ £ ) 



»7 (r, k 2 )f] b (T, k 2 ) - f) a (T, k 2 )rj h {T, k 2 



-W . 



(33) 



as can be checked by an explicit calculation. The right hand side of the previous 
equation is nothing but the Wronskian W of the pair of solutions r\ a ,b and is 
independent of r in the case of eq. (29). Let us denote by W the value of the 
Wronskian for the pair of solutions rj a ^. The discontinuity of d T 'H(r, r', k 2 ) 
across t' = r is therefore equal to W, which means that the second time deriva- 
tive indeed contains a term W 5(t — t'). Finally, from the obvious property 

Q(t, r, k 2 ) = , (34) 

one easily sees that H(t,t' ,k 2 ) satisfies the boundary condition of eq. (28). 
Therefore, W 1 H(t, t', k 2 ) is the propagator we are looking for: 

g(t j h) _ j¥A max(r, r% fc 2 )n(min(r, r% -(3/2, k 2 ) 

1 ' ' ' wn(i3/2,-p/2,k 2 ) ■ 1 ' 

In general, the solutions r\ a ,b will not be found analytically for a non-zero k, and 
will have to be found numerically. 



3.4 Calculation of the functional determinant 

As we have already explained, we need to calculate the determinant that appears 
in eq. (18), det (A*[(/> c ]), to order zero in the fluctuation £(x) of the boundary, 
i.e., det (A* [(/>o]) . Integrating by parts the kinetic term, the Euclidean action 
can be rewritten as : 



S E [</>o + rf\ 



(d 4 x) 



^rjD^+^m 2 ^ + U"(cj)o)v 2 



(36) 



Notice that the integration by parts does not introduce any boundary term here, 
thanks to the boundary condition obeyed by rj (see eq. (7)). Therefore, we have 
for the operator A* the following expression : 

A*[4>o\ t ,x;t' ,y = 



84>{t,x)84>{t' ,y) 



4>=<t>o 



= S(T-T')8(x-y)[n E +m 2 + U"(Mr))] ■ (37) 
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Notice that here we have already written <j>o explicitly as a field that depends 
only on time (because we are calculating the determinant only at order zero in 
the fluctuations of the boundary). Thus, we can perform a Fourier transform 
with respect to space, and use fc instead of x. The eigenvalues 9 i and eigen- 
functions r)i of the operator A* [<j>o] are functions t](t, x) that obey the following 
system of equations 6 : 

d 2 -(m 2 + k 2 )- U" (Mr)) ] Vi(T, k 2 ) = 9i r,i(r, k 2 ) , 
Vas, Vl (-p/2,k 2 ) ^ Vl {(3/2,k 2 ) -0. (38) 

This equation is of the same type as eq. (29), the only difference being that k 2 
is now replaced by k 2 + gi. Therefore, it has two independent solutions that are 
given by r] a (T, k 2 + gi) and ^(t, k 2 + gi), and its general solution can be written 

as: 

77 4 (r, fc 2 ) = C a j] a (r, k 2 + gi ) + C b r] b (T, k 2 + 9i ) , (39) 

where C a ,b are two integration constants. In order to have a non-zero r\i that 
obeys the required boundary conditions, we need to have the following property: 

r] a (-P/2, k 2 + gM(3/2, k 2 + 9l ) = ri a (0/2, k 2 + gM~p/2, k 2 + 9i ) . (40) 

This equation determines the allowed eigenvalues gi. This equation can also be 
written as : 

O(/?/2,-/3/2,fc 2 +. 9l ) = 0, (41) 

where fi has been introduced in eq. (31). The determinant of the operator A* 
is of course obtained as the product of its eigenvalues: 

detA*[^»o,fc 2 ]= I] 9- (42) 

g\n(p/2 ,-/3/2,fc 2 +9)=0 

(We denote by A*[(j)Q,k 2 ] the restriction of the operator v4*[</>o] to field fluctu- 
ations of Fourier mode fc.) If we denote by z n the (possibly complex) zeros of 
the function Sl(0/2,-/3/2;z), then the solutions of Q(/3/2, -(3/2, k 2 + g) = 
are the numbers g = z n — k 2 . Therefore, we can write 

dot A*[Mk 2 } =\{{z n -k 2 ) , (43) 

n 

where multiple zeros are repeated as many times as needed in the product. 
The right hand side of this equation is an entire function of fc 2 , that obviously 
vanishes at all the z n 's. Since Q(fi/2, —j3/2, fc 2 ) shares the same property, there 
exists an entire function p(k 2 ) such that [37] : 

dot A* [</.„, fc 2 ] = ft(/?/2,-/3/2,fc 2 ) e p(fe2) . (44) 



"For many more informations about properties of Hill's equations and their solutions, the 
reader may consult [36]. 
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From eq. (31) we see that, if we chose the functions rj a and rjb in such a way 
that their value at r = —(3/2 is independent of k 2 , then the limit 

Jim \n({3/2,-p/2,k 2 )\e- M ^ (45) 

is bounded for any M > (3. By Hadamard's theorem [37] , we conclude that the 
function p(k 2 ) is a constant 7 . The constant factor e p must in fact be propor- 
tional to the inverse of the Wronskian of the pair of solutions r\ a and rjb that 
we are using, e p = const//3VF, because the determinant must be independent 
of this choice (the factor f3 in the denominator has been included by hand in 
order to have a dimensionless determinant). The normalization constant can be 
absorbed as an overall multiplicative constant in Z. 

Finally, the determinant of A* [fa] is obtained by multiplying the previous 
result for all fe's, which gives : 

. ,*r, 1 j t f d3k , fflW/2,~/3/2,k 2 )\ 

det A* [ M = exp V J — In ( ^ ) ■ (46) 

In order to see how the volume V appears in this formula, it is useful to con- 
sider first that the system is in a finite box, and to rewrite the sum over the 
corresponding discrete Fourier modes as an integral. 



4 Integration over the boundary fluctuations 

4.1 Expansion of the classical action 

The final step in the analytic part of this calculation is to calculate the functional 
integral over the fluctuation £(x) of the boundary in eq. (18). Before doing this 
integration, we must expand the classical action S E [<j> c ] to quadratic order in £, 
using the expansion of eq. (17) for <fi c . We have: 

S E [d> c ] = S E [<M + S^S E + S^S E + O(e) ■ (47) 

Notice that we could be in trouble because a priori we must keep (j>2 - the term 
of order £ 2 in the classical solution <j> c - in the second term of the right hand 
side, which would be much more difficult to obtain. We will not need this term 
however, because 0o is an exact solution of the classical equations of motion. 
Indeed, one can write 

S^S E = J (^x), \ [-Ur)+m 2 Mr) + U'(Mr))\ {4>i{x) + <h(?)) 



d A x 



Mr){M*) + h (a;)) • (48) 

— p/2 



7 Strictly speaking, this result only proves the independence of the function p with respect 
to fc 2 , but it does not exclude a dependence on the other parameters of the problem: the 
mass m and the coupling constants contained in the potential U(<j>). However, as is clear 
from the operator whose determinant we are calculating, this dependence only arises from the 
combination (m 2 + U"((/>o(t)) which means that it can enter in the final result only via the 
solutions rja and rj b , i.e., via the function Q. Thus, the prcfactor cxp(p(fc 2 )) cannot contain 
any implicit dependence on these parameters. 
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The integrand in the first term of the right hand side vanishes identically because 
4>o(t) obeys the classical equation of motion associated to the action S E . The 
second term - a boundary term - can be rewritten as follows : 



/ 



d A x 



Mt)(M*) + fate)) 



t=+P/2 
T=-f3/2 



0o(+/J/2)-0o(-0/2)J J d 3 xti(x) , 

(49) 

and it vanishes because the fluctuation of the field at the boundary has a 
vanishing average. The second order variation of the action - the third term in 
the right hand side of eq. (47) - can be written as 

* (2) S E = \J (d 4 x) E [(a^ 1 (x))(a^ 1 (x))+m 2 ^(x) + L/"(0o(r))^(x)) 



(d 4 x) E a M M*WMx) 



+ \ J (d 4 x) E M*) p B +m 2 + U"(Mr))} M*) ■ 



(50) 



The integrand of the second term vanishes because of the equation of motion 
obeyed by the field <j>i(x). Therefore, the second order variation of the classical 
action comes entirely from the boundary term 



5&S E = \j d 3 x[Mr,x)d T Mr,x)} T T Zt%l ■ 



(51) 



By rewriting this integral in momentum space, and by making use of the bound- 
ary condition obeyed by </>i(t, x) and of eq. (26), we can rewrite this as follows 8 : 



(52) 



where we denote 



h8/2- 

-P/2. 



r=+P/2 
t=-0/2 



(53) 



C(k)= a r a T /G(r,r',fc) 
Therefore, the Gaussian functional integral over £ leads to the following result: 



nw) 

k 



e - S e [M y/PC(0) CX P 



V 



d 3 k 



■In {pC{k)) 



(54) 



The factor y/ (3C(0) in the numerator serves to remove the contribution of the 
zero- modes, i.e., the functions £(x) that are constant, since these are taken into 
account in the quantity ipo. Factors of j3 have been introduced in order to make 
the arguments of the log and of the square root dimensionless. 



8 In order to obtain this formula, we use the relation 

=+(3/2 _ n(T,-§,k 2 ) + n(§,T,k 2 ) 



[d T ,G(r,r',k)] T T ,: 



-0/2 



fXf,-§,fc 2 ) 



Therefore, this quantity is equal to 1 at r = ±/?/2. 
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4.2 Calculation of C(k) 

The quantity C(fe) defined in cq. (53) involves the calculation of two derivatives 
of the Green's function evaluated at the boundaries. This may pose a problem 
because the derivative of G is not continuous at coincident points. It is crucial 
to note that eq. (53) imposes a very definite order when taking the limits r, r' — ► 
±/3/2. This leads to an unambiguous expression for C(k) : 

C(k) = lim d T d T/ G(T,T , ,k)+ lim <9 t <9 t ,G(t, t', k) 



(9/2 



- lim d T d T ,G(T,T',k) - lim <9 T <9 r /G(T, r' ,k) . (55) 

T '^ + /3/2 T l-,-p/2 

From eq. (35), we see that, depending on the order of r and r', the double 
derivative of the propagator reads : 



d T d T ,G{T, T ' ,k) = - 
a r a r ,G(r,T',fc) = - 



d T n(/3/2, t, fc 2 )d T ^(r', -/?A fc 2 ) 
Wfi(/?/2,-/3/2,fc 2 ) 

d T >tl{0/2, t', k 2 )d T il( T , -(3/2, fc 2 ) 
W n(/3/2,-(3/2,k 2 ) 



if t' < t , 



if r' > r . 



(56) 



Using the explicit form of fi(r, t', fc 2 ) given in eq. (31), a straightforward calcu- 
lation gives : 

'Ar? a (fc 2 ) A?y a (fc 2 
y A Vb (k 2 ) Ar) b (fc 2 



det 



C(k) = 



det 



7? a (§,fc 2 ) r? a (-f,fc 2 ) 
r? 6 (f,fc 2 ) Vb (-lk 2 ) 



(57) 



where we denote 



A7? a ,t(fc 2 ) = r) a , b (T,k 2 ) 



T=+f3/2 
r=-/9/2 



A?) a ,h(fc 2 ) = fla >b (T,k 2 ) 



r=+P/2 

T=-0/2 



(58) 



Notice that the form of C(fc) given in eq. (57) makes obvious the fact that C(fc) 
does not depend upon the choice of the two solutions rj a and rjb that one takes, 
as long as they are linearly independent. Indeed, the coefficients C(k) are a 
property of the classical action itself, and should be independent on the basis 
chosen for the fluctuations around the classical field. 
If we take two solutions rj a and rj b such that 



ry a (-/?/2;fc 2 ) = l, 7? a (-/?/2; fc 2 ) = , 
%(-/3/2;fc 2 )=0, 7/ b (-/3/2;fc 2 ) = l//3 



then 



C{k) 



r? a (/3/2)+/?^Q3/2)-2 
W/2) 



(59) 
(60) 
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We will suppose that (m 2 + U") is positive 9 . In this case, one can easily show 
from (29) that rj a and rj b are monotonically increasing in [—(3/2, /3/2]. This 
implies that C(k) > and S^S E > 0, which means that the fluctuations of 
the boundary field always increase the value of the action compared to the 
configuration with a uniform boundary condition. This can be seen as an a 
posteriori justification for the choice of expanding around configurations with a 
uniform boundary condition; indeed, such configurations have a smaller action 
than those with fluctuations of the boundary condition, and thus are the leading 
contribution to the partition function. 

4.3 Diagrammatic interpretation 

The Gaussian integration of exp(— S E [</> c ]) over the fluctuations of the field on 
the time boundary also corresponds to some one loop corrections. To begin 
with, let us recall the obvious fact that the classical action S E {<j> c ] only contains 
terms that are quadratic or quartic in the classical field <j) c . Moreover, we have 
already seen at the end of section 2 that the classical field <p c is the sum of all 
the tree diagrams with one external leg, terminated on the other side by the 
boundary field ip (sec figure 2). Thus, S E [(j) c ] is a sum of tree diagrams that 
have no external legs, with the boundary field <p at the endpoints of the tree. A 
typical diagram of that sort has been represented in figure 3. 




Figure 3: Diagrammatic expansion of the classical action S E [</> c ] in terms of the 
boundary value of the field (black dots). 

At this point, these diagrams represent the classical action for an arbitrary 
field ip as the boundary condition. Writing <p(x) = p a + £(x) and doing a 
Gaussian approximation means that, for each diagram like the one displayed 
in figure 3, all the black dots except two of them are replaced by a uniform 
boundary field ipo and the remaining two are replaced by the fluctuation £(x) of 
the boundary. Then, integrating out the field £ means that the endpoints where 
the £'s are attached are linked together, thereby forming a loop. To this loop 
can be attached an arbitrary number of tree diagrams terminated by ipo : each 

9 In other terms, the spectrum of the semiclassical propagator has no bound states [36]. 
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of these trees is a contribution to <j>o{r), the classical solution with boundary 
value (pa. 

Thus, we conclude that the terms resulting from the Gaussian average over 
the fluctuations of the boundary field are also 1-loop contributions in a back- 
ground made of the field </>o (r) . These terms are therefore on the same footing 
as the terms included via the determinant det (A* [0o]). Moreover, this analysis 
of the diagrammatic content of our approximate expressions confirms the self- 
consistency of these approximations: it would have been inconsistent to keep 
Gaussian fluctuations of the boundary in det (A* [0 C ]), because by doing this we 
would include two- loop terms in the background field 4>q. 

As we shall see in section 5, another consistency check of our final formula can 
be made based on the structure of its ultraviolet divergences: it contains exactly 
the divergences one expects of the 1-loop effective action in the background field 
4>o (r), and is thus straightforward to renormalizc. It is important to realize that 
we need both the 1-loop corrections coming from det (A* [<p c ] ) , and those coming 
from the Gaussian integration over the fluctuations of the boundary field in 
order to reproduce the usual pattern of 1-loop ultra-violet divergences. Failing 
to include one of the types of terms, one would have spurious divergences that 
could not be removed by the usual renormalization procedure. 



4.4 Final formula for the partition function 



Collecting everything together, we can write the following formula for the (non- 
renormalized) partition function: 



+00 

/ 



dipo e 



-S E [<j> ] 



#(0) 



x cxp — - 



V f d 3 k 



(2tt) 3 



hi 



1 

W 



Ar, a (k 2 ) Ar, a (k 2 ) 
A Vb (k 2 ) Af] b (k 2 ) 



(61) 



which is valid for arbitrary choices of fc 2 -independent initial conditions. Indeed, 
the ratio of the determinant and the Wronskian inside the logarithm does not 
depend on any particular choice for the two solutions r\ a and 77;,. In practice, 
one can take advantage of this freedom in order to simplify the numerical cal- 
culations. In particular, for the initial conditions defined in (59) we have 10 



1 

W 



A Va (k 2 ) 
A Vb (k 2 ) 



Af, a (k 2 ) 
Ar, b (k 2 ) 



= ry a (/3/2,fc 2 )+/3^(/3/2,fc 2 )-2. 



(62) 



Thus, we have obtained a fairly compact formula that resums (in the Gaussian 
approximation) the fluctuations around the classical solution and the fluctua- 
tions of the boundary condition. At this stage, the calculation only involves 



10 This can be further simplified when the potential U{(j>) is an even function of 1 
case, we have rj a (fi/2, k 2 ) = /3rj b (/3/2, k 2 ). 



In this 
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solutions of some ordinary differential equations, which is in principle straight- 
forward to obtain numerically. For each tpo, one must determine the following 
quantities: 

1. the classical solution <po{r), 

2. the classical action S E [<j>o], 

3. for each k 2 , two independent solutions r) a (r, k) and r)b(r; k) of the equation 
of fluctuations around the classical solution (3 (t). 

Notice that all the quantities that depend on k in fact only depend on \k\. This 
means that the integration over k is in fact a one dimensional integral. 



5 Renormalization 

Our final expression, eq. (61), is plagued by ultraviolet divergences if taken at 
face value. These divergences arise from the integration over the momentum k 
in the second line. It is in fact easy to convince oneself that these divergences 
can be dealt with by the usual 1-loop renormalization procedure. In order to 
see this, one must write the solutions r] a and rjb as series in the interaction term 
U"((f>o) with the background field. Indeed, if we denote by rf™} the term in rj a ^ 
that has n powers of U"((j>o), we have the following relations : 

(d 2 T = , 

From these equations, one can see that ^ b +1 ^ has an extra power of 1/fc 2 at 

large k compared to rfc 1 }. Thus, we expect that only a finite number of terms in 
this expansion will actually contain ultraviolet divergences. To check this, let us 
calculate explicitly the first three terms in the expansion of the right hand side 
of eq. (62). The solutions rf®\ that obey the boundary conditions of eq. (59) are 
given by : 

^°)(r,fe 2 )=cosh(c fc (r+^)) , 

_(0) (T)fc2) = sinhM, + f)) _ (g4) 

Notice that these Oth-order solutions already saturate the boundary conditions 
at r = —13/2 in eq. (59). Thus, the higher order terms in rj a b should vanish and 
have a vanishing first time derivative at r = —(3/2. In order to find these terms, 
it is useful to first construct a Green's function G°(t, r',fe 2 ) of the operator 
d 2 — c<j| that obeys the following conditions : 

(d 2 -cj 2 k )G^(r,r',k 2 ) = S(r-r'), 

G^(r = -^,r',fe 2 ) = , d T Gt(T = -^,T>,k 2 ) = 0. (65) 
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It is straightforward to check that the propagator obeying these conditions is 
given by 



which is nothing but the retarded Green's function of d 2 — uj\. With this Green's 
function, one can write 



r+0/2 

dr'G^T 1 r'^)U''(Mr'))v { a n l(r'^). (67) 

J -fill 



Notice that, since the classical solution 4>q(t) does not depend on space, the 

relationship between rf^t^ an< ^ ^i™fc 1S l° ca l i n k. 

At this point, it is a straightforward matter of algebra to obtain rj a b up to 
second order in U" . We obtain 



. f +f3/2 dT , U"(Mr')) { 1 

J-fi/2 2wfc 2 

1 f +P/2 , e -2o> fc |r'-r"| 

2 



+/3/2 

■fi/2 2cj fc 



0/2 dT ' dT " (2c k y U "W» u "(Mr")) 



+0(e-^) + 0((U") 3 )j . (68) 

Inside the curly brackets, we have dropped all the terms that would go to zero 
exponentially when \k\ — > +oo. Indeed, these terms do not contribute to the 
ultra-violet divergences we are studying in this section. In this expression, we 
recognize the time-ordered propagator, which reads 

„ p-Wfc|T-T'| 

G°(r,r',fc 2 ) = ^ . (69) 

It is a remarkable feature of eq. (62) that, while having a fairly natural expression 
in terms of a retarded propagator, it can be rearranged as an expression involving 
the time-ordered propagator (at least for the terms that will contribute to the 
ultra-violet divergences) . 

The terms that appear in the curly bracket in eq. (68) have a fairly simple 
interpretation in terms of Feynman diagrams. For a scalar theory with a </> 4 
coupling, the first non-trivial term can be represented as 




U"(Mr'» _ y ; ( . 0) 
fi/2 2w fc 
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Notice that, in this expression, l/2u>k is the equal-time value of the time-ordered 
propagator. Similarly, the term on the third line can be represented as 

1 f +l3/2 e -2o> fc |r'-r"| \ / 

-d,„ wi ^r r * (/ » B " w ' </ » - XZX- 

(71) 

The second term on the second line of eq. (68) would be represented by a graph 
made of two disconnected components, each of which is given in eq. (70) (the 
factor 1/2 is the symmetry factor that results from the possibility of exchang- 
ing the two connected components). In fact, when we take the logarithm (as 
required by eq. 61), these disconnected contributions simply drop out : 

l ln r. ,0 ,2, , ^ ,0 ,2, J _ fa , 1 f +P/2 ^ U"(Mr')) 



2 



V a {^) + ^/-^)-2]=?f + \ J dr 



2u k 



1 /-+/3/2 p -2w k \r'-T"\ 

-s J_ m dT ' dT " u " {Mt ' )] u " (Mt " )] + (72) 

One can check that the cancellation of the disconnected terms when one takes 
the logarithm is in fact quite general, and works to all orders. Finally, when we 
integrate over k, the first term gives the usual zero-point energy, and the next 
two terms are the first two non-trivial terms of the zero temperature 11 1-loop 
effective action (for this, it was important to be able to rewrite the expression in 
terms of time-ordered propagators). All these terms are ultra-violet divergent. 
If calculated with a momentum cutoff A, they behave respectively as A 4 , A 2 , 
and ln(A), if there are 3 spatial dimensions. All the higher order terms in the 
expansion in powers of U" are ultraviolet finite, because they have at least one 
extra power of 1/fc 2 when \k\ — » +oo. 

This identification tells us that, in order to renormalize our final expression, 
we must follow the following procedure : 

1. subtract the "zero point energy" in ln(Z), i.e., subtract /3wfc/2 from the 
integrand in the integration over k, 

2. add the one-loop counterterms to the classical action S E [(pa], and simul- 
taneously regularize the integration over k. 

Notice that the regularization scheme employed for calculating the counterterms 
must be identical to that used when computing the integral over k. Thus, a 
regularization by an ultra-violet cutoff seems the most convenient method here. 
Once the above two steps have been carried out, one will have a A dependent 
expression that tends to a finite result when A — ► +oo. 

This expression of Z is free of any ultra-violet divergence. But, naturally, it 
is now expressed in terms of couplings and masses that are scheme dependent 
(because one must chose a particular renormalization scheme 12 in order to define 



n We recover the well known fact that, if a theory is renormalizable at T = 0, it is also 
renormalizable at finite T, with the counterterms evaluated at T = 0. 

12 The renormalization scheme should not be confused with the regularization scheme. 
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uniquely the counterterms that are added to the classical action) . The standard 
procedure at this point is to express other physical quantities in terms of the 
same scheme-dependent parameters, and to eliminate them in order to have 
relationships that involve only physical quantities. 

6 Conclusions 

We have derived a semiclassical approximation for the partition function of a 
system of scalar fields in the presence of an arbitrary single-well interaction po- 
tential. In the path-integral formalism, the partition function is an integral over 
periodic configurations in imaginary time, and is dominated by classical tra- 
jectories. The non-perturbative information contained in the classical solutions 
serves as the starting point for this semiclassical approximation. 

Euclidean classical solutions are usually not known for arbitrary (periodic) 
boundary conditions. However, by first considering classical solutions that cor- 
respond to a spatially independent boundary condition (finding these special 
solutions amounts to solving an ordinary differential equation), one can con- 
struct approximate classical solutions obeying arbitrary boundary conditions in 
a systematic fashion. We have calculated the contribution of quantum fluctu- 
ations around those classical solutions in a self-consistent scheme. Our final 
formula for Z admits a simple expression in terms of two independent solutions 
of the equation of small fluctuations around the classical solutions, and is thus 
easily amenable to a numerical evaluation. Despite its simplicity, our expres- 
sion treats exactly the mean value of the field on the boundary, no matter how 
large. Moreover, we have shown that this expression is renormalizable by the 
subtraction of the standard one-loop counterterms, and by the subtraction of 
the free-field energy. 

The formula we have obtained for the partition function is non-perturbative 
in the sense that it resums the interactions to all orders for the configurations 
where the mean value of the field on the boundary is large. This can be seen 
by investigating which classes of diagrams of the usual perturbation theory are 
taken into account in our approach. We expect that thermodynamical prop- 
erties derived from this semiclassical expression for Z will be valid in a wider 
domain in the parameter space (T, {A}) (where {A} represents the coupling con- 
stants) as compared to results obtained from the plain perturbative expansion. 
We are currently investigating in detail the case of a theory with a \<p 4 cou- 
pling. Results, including a detailed comparison with those obtained by other 
resummation schemes, will be presented in a future publication. 

Natural candidates for a direct application of the result derived in this paper 
are condensed matter systems containing scalar order parameters, such as den- 
sity or magnetization. Extensions to potentials with more that one minimum, 
and other field theories can also be pursued. 
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A The free case 



/•/J/2 

S[</>] = / d 3 xdr 



2 ^<W + -m 



The action in the free theory is given by 

"1 

d"xdr 

-13/2 

leading to the equation of motion: 

[d 2 T - m 2 ]0 o = . 
The classical solution satisfying <j>o(—(3/2) = (f> ((3/2) 



2-2 



ip IS 



cosh m(r + 13/2)) + (1 .^j^ sinh(m(r + (3/2)) 



(73) 



(74) 



(75) 



It is easy to show that S^^o] = onp\, with a = mV (cosh(f3m) — l)/sinh(/3m), 
where V is the volume. Following our main result, we need two solutions of 

d 2 - (m 2 + k 2 )]r) = , (76) 

obeying eq. (59). We have already seen these solutions in eq. (64). We obtain 

rj a (/3/2, k) + (37j b (/3/2, k)-2 = 2(cosh(/? Wfe ) - 1) (77) 

= (1 - cxp(-/?w fe )) 2 exp(/3w fe ) 

and 

2m(3 (cosh(/3m) - 1) 2a(3 



f3C(0) 



Finally, we have 



2a/3 
~V~ 

2^(3 



■ exp 



exp 
d 3 k 



sinh(/3m) 



/ ^ e ~ avi exp [- y J H 



V 



) + 



(3uJk 



2 J 



(78) 



(79) 



(80) 



that is (up to an overall constant) the known result for the harmonic oscillator 
(not yet renormalizcd). We see that our approximation scheme leads to the 
exact result in the case of the free theory. Naturally, this is due to the fact 
that, in the absence of any interactions, the Gaussian approximation represents 
exactly the fluctuations in the system. 
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